
#include "PyrolysisHeatDiffusionKernel.h"

registerMooseObject("TrilobitaApp", PyrolysisHeatDiffusionKernel);
template<>
InputParameters validParams<PyrolysisHeatDiffusionKernel>()
{
  InputParameters params = validParams<Kernel>();
  params.addClassDescription("Compute heat diffusion kernel ");
  return params;
}

PyrolysisHeatDiffusionKernel::PyrolysisHeatDiffusionKernel(const InputParameters & parameters) :
  Kernel(parameters),
  _ks(getMaterialProperty<Real>("ks"))
{

}

Real PyrolysisHeatDiffusionKernel::computeQpResidual()
{
//	std::cout << _k[_qp] <<std::endl;

  Real r(0);
  r = _ks[_qp]*_grad_u[_qp] * _grad_test[_i][_qp];
  return r;
}

Real PyrolysisHeatDiffusionKernel::computeQpJacobian()
{
  Real jac(0);
  jac = _ks[_qp]*_grad_phi[_j][_qp] * _grad_test[_i][_qp];
  return jac;
}
